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Abstract 

A physically defined effective charge can incorporate quark masses analytically at the flavor thresholds. 
Therefore, no matching conditions are required for the evolution of the strong coupling constant through these 
thresholds. In this paper, we calculate the massive fermionic corrections to the heavy quark potential through 
two loops. The calculation uses a mixed approach of analytical, computer-algebraic and numerical tools includ- 
ing Monte Carlo integration of finite terms. Strong consistency checks are performed by ensuring the proper 
cancellation of all non-local divergences by the appropriate counterterms and by comparing with the massless 
limit. The size of the effect for the (gauge invariant) fermionic part of ay(q 2 ,m 2 ) relative to the massless case 
at the charm and bottom flavor thresholds is found to be of order 33%. 
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1 Introduction 



In analogy to Quantum Electrodynamics, the heavy quark potential has been of interest in QCD from very early 
on [16, 34, 35, 41, || |lj as a model for the physical definition of the strong coupling constant [29]. Since it 



represents a potentially measurable quantity and gives naturally rise to a physical effective charge ay |29J, it is 



very interesting to study the QCD flavor thresholds in such a system [30] as the fermionic corrections are separately 
gauge invariant. 



In the MS and the MS schemes, the running of the coupling constant, by construction, does not know about 
masses of quarks and since the couplings are non-physical, the Appelquist-Carazzone |33j decoupling theorem is 
not applicable. One has to turn to effective descriptions which match theories with m massless flavors onto a 



theory with m — 1 massless and one massive flavor at the "heavy" quark threshold 15]. In this way, the 

dependence on the dimensional regularization mass parameter [i is reduced to next to leading order effects by 
giving up the analyticity of the coupling at the flavor threshold [^, 39, 37. 38, [14], |9[. 



While this procedure of matching conditions and effective descriptions is certainly workable, from a theoretical 
standpoint it would be advantageous to have a physical coupling constant definition which is analytic at thresholds. 
In addition, as a physical observable, the total derivative with respect to the renormalization scale [i vanishes. 
Such a system is given by identifying the ground state energy of the vacuum expectation value of the Wilson loop 



as the potential V between a static quark-antiquark pair in a color singlet state [16, 41, 21 1: 



V(r, m 2 ) = - hm hog(0\Tr {P exp (j dx^T^ }|0) (1) 

where r denotes the relative distance between the heavy quarks, m the mass of "light" quarks contributing 
through loop effects and T a the generators of the gauge group. It is then convenient to define the effective charge 
ay (q 2 , m 2 ) as 

2 a An^ay^rrfl 
q 2 

in momentum space. The factor Cp is the value of the Casimir operator T a T a in the fundamental represen- 
tation of the external sources and factors out to all orders in perturbation theory. As one is free to choose the 
representation of the external particles, we obtain the static gluino potential by adopting the adjoint representation. 



The massless case was recently calculated in Rcf. [22] and in this paper, we will give all the two loop fermionic 



contributions to ay (q 2 , m?) for all perturbative values of the momentum transfer q 2 = q$ — q 2 = — q 2 > and for 
arbitrary values of the fermion mass m. In this context we are only interested in the two loop contributions to the 
potential in the effective Schroedinger equation for the heavy particles. This implies, for instance, that not always 
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the whole diagram contributes to the potential as certain parts can already be reproduced by the exponentiation 
of lower order diagrams. The necessity for this subtlety has its origin in the exponential present in Eq. |l[ For a 
detailed discussion, see Ref. [41]. 



It is also important to note that the results of massive two loop integrals presented in this work are also relevant 
for the related problem of quark threshold production. For this application, though, it would be necessary to treat 
also the occurring imaginary parts of the integrals numerically as pole terms will contribute for timelike momentum 
transfers at the production threshold q 2 = 4m 2 . A promising approach for this treatment might be the recently 
suggested Taylor expansion of integrands around threshold |lg] by determining large and small scales in the 
problem. The heavy quark approximation eliminates the possibility of timelike momentum transfers in this work 
so that we do not need to worry about pole terms numerically. Nevertheless, we also list the contributions needed 
in this case for all integrals. 

The paper is outlined as follows: 

In section [2] we list all the occurring two loop contributions explicitly in the Feynman gauge and with the 
usage of heavy quark effective Feynman rules for the external sources. In section ^ the unrenormalized results 
for the two loop corrections are given in terms of two loop scalar integrals, for which explicit expressions are 
listed in appendix [B]. Section || contains all the required counterterms in the M£-renormalization scheme and it 
is shown that all non-local divergences cancel. The renormalization constants obtained are given explicitly and 
checked with the known results. Section || contains numerical results which demonstrate that the massless limit is 
obtained correctly and display the effect of including the mass terms for the charm and bottom flavor thresholds. 
In section || we make concluding remarks and indicate future lines of work with the presented results. Appendix 
|A| , finally, lists all the reductions from tensor to scalar integrals needed for the results displayed in section ||. 

2 The Two Loop Corrections 

In this section we present the non-Abelian contributions to the heavy quark potential that constitute the new 
results of this work. They are depicted in Fig. The QED like diagrams, which need to be modified by their 
respective color factors, have been known for a long time || and can also be found in Refs. [14, 0, 0, 27] for instance. 



They are given here as well because we would like to be able to separate non-Abelian and Abelian contributions 
to the potential. It has been observed before []l4| that their respective threshold behavior can be quite different. 
These diagrams, together with effectively "one loop" diagrams are given in Fig. |3|. The weighted sum of all the 
graphs shown, modulo terms already generated by the exponentiation of the lower order Born and the one loop 
vacuum polarization diagram, give the complete gauge invariant fermionic corrections to the heavy quark potential 
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Figure 1: The Feynman rules for heavy quark effective theory used in this work for the source propagator and the 
source gluon vertex. For anti sources one has to make the replacement v — * —v. The is prescription is the same 
as for the usual fermion propagator. 

at two loops in the Feynman gauge. The choice of this gauge simplifies the calculation because the decomposition 
into scalar two loop integrals is easier and it also reduces the three gluon vertex correction graph to zero in the 
heavy quark effective theory. Below we list all contributions at the two loop level. The abbreviations stand for 
gse = gluon self energy, vc = vertex correction, cl = crossed ladder and olvc = one loop vertex correction. 
In the heavy quark limit we use the source gluon vertex and source propagator Feynman rules of heavy quark 
effective theory [19, fn which are given in Fig. ffl. 



With these, and taking v^ = (1, 0, 0, 0) and qo = for the purely spacelike momentum transfer q, the two loop 
diagrams of Figs. |2] and || read in the Feynman gauge (summed over the external color degrees of freedom and 
including a symmetry factor of ^ for the first three amplitudes): 
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(9) 
(10) 

(11) 
(12) 



It should be noted that in our case there is no need for an is prescription in the denominators of Eqs. || through 
|9| as the spacelike nature of the physical momentum transfer only leads to purely real integrals and no unambiguous 
pole terms occur in the denominators of those diagrams. This feature also simplifies the Monte Carlo integration 



of the finite parts of the contributing graphs. The three graphs 10, 11 and 12 display infra-red divergences which 
cancel in the sum. The one loop vertex correction graph M i vc vanishes in dimensional regularization, however, 
is needed to ensure the proper cancellation of infra-red divergences. 

The color factors given are not always the full color factors. Only those contributing to the potential are listed. 
The Casimir invariants [^] for a general SU(N) group are defined by 

N 2 - 1 

Ca = N , C F ^^- (13) 



Furthermore, Tr{T a T b } = T f 5 a ' b = \5 a > b . The color factor for M VCl includes the sum of the graph shown in 
Fig. |2] plus the term stemming from the fermion momenta reversed contribution. Only the sum is proportional 
to Ca, the other terms vanish according to Furry's theorem, as is the case in QED. For QCD, the crossed ladder 
diagrams do contribute as they contain a color factor proportional to Cp — Cp ^ A , whereas the straight ladder 
graph has a color factor proportional to Cp only. This will be expounded on in section |3.1| . In QED, the sum 
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gse^ gse 2 : gse 3 : 




Figure 2: The non-Abelian Feynman diagrams contributing to the massive fermionic corrections to the heavy 
quark potential at the two loop level. The first two rows contain diagrams with a typical non-Abelian topology. 
Double lines denote the heavy quarks, single lines the "light" quarks. Color and Lorentz indices are suppressed in 
the first graph. The notation for the remaining digrams is analogous. The last line includes the infra-red divergent 
"Abelian" Feynman diagrams. While the topology of these three diagrams is the same as in QED, they contribute 
to the potential only in the non-Abelian theory due to color factors CfCa- In addition, although each diagram 
is infra-red divergent, their sum is infra-red finite. 
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gse 4 : gse 5 : 2vp: 




Figure 3: The infra-red finite Feynman diagrams with an Abelian topology (upper line) contributing to the 
massive fermionic corrections to the heavy quark potential at the two loop level plus diagrams consisting of one 
loop insertions with non- Abelian terms (lower line). 

of all vertex, ladder and crossed ladder Feynman diagrams are equivalent to the iteration of the potential in the 
Schroedinger theory. av QED and the effective coupling || differ, therefore, only at three loops due to light by light 
scattering contributions. 

3 Unrenormalized Results 

The two loop integrals needed for the expressions of Eqs. | through | are treated in separate ways in this work 
depending on whether or not they contain two or more internal fermion lines. In the former case we integrate 
the fermion loop first as will be explained below. For the vertex correction contribution Ai VCl we integrate the 
fermion loop analytically as well with all the Lorentz indices projected to zero and then proceed with additional 
Feynman parameters for the remaining loop integration. 

The two point functions A4 gsei , A4 gse4 and M gse5 are treated in a completely different manner as the above 
techniques would now be too cumbersome. We project the complicated tensor structure onto scalar quantities 
as described below and then proceed with an algebraic reduction into scalar two loop integrals. This reduction 
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is programmed in FORM [13| and details are presented in appendix |A[ The resulting scalar integrals are then 
evaluated by employing standard Feynman parameter techniques and explicit results are listed in appendix [B|. 
Overall results for the various amplitudes are obtained by expanding the n-dimensional results around e = with 
MAPLE. It is important to notice, given the complexity of the calculation, that the translation into FORTRAN 
code was also performed by MAPLE, thus dramatically reducing the chance of accidental mistakes. The evaluation 
of finite parts is done with the Monte Carlo integrator VEGAS [25[. 

For the two point functions we use the following decomposition into transverse (t) and longitudinal (I) com- 
ponents: 



n^(.)-(^-^)n t (, 2 ) + ^n,(, 2 ) 



from which it follows that in n = 4 — e dimensions 



(14) 



n — 1 



<zV n ( x 



(15) 
(16) 



With this notation and the heavy quark effective Feynman rules depicted in Fig. |l] we arrive at 



M gsei = (q) ,i = {1...5} 



(17) 



The result of the decomposition for the transverse component of the gluon self energy graph A4 gsei , using the 
relations given in appendix [A|, reads 



n, 1 



ig 4 C A T F 
4(n - 1) 



8 20 



n- 



+ 16 ( q 2 - q 2 ^ - m 2 ) T 1235 + q 2 (4n - 6) T 23 45 + q 2 (2n - 4) B 12 B 45 - 8q 2 m 2 T 12U5 + 8A 2 B i5 



) (a 2 B 12 , - m 2 T 12 , 35 ) + (An - 10) T 235 + (8 - 4n) A 2 B 12 + u| - ^ T 135 



+ 



m" I / 8 20\ ( 9m m „ „ 1 ,-, 



71- 



"7 Ti 2 / 35 + T 2 / 35 - A 2j B 12 / H g A 2 + 16T 235 



77? '■ 



77- 



28\ 1 1 



-n^77i 2 Ti 2 / 3 5 + nT 235 - re^Ti 35 - 4nA 2 B i5 + n^A 2 B 12/ - 8m 2 T 2U5 



-4m 2 T 2355 + 4T 135 + g 2 (4777 2 r 2345 5 + (4 - 2n)^ 2 C 4 55 - nT 2Mb ) + ^ (4m 2 (T 135 - T 235 ) 
2 

+n- (m 4 Ti 2 / 3 5 - m 2 Ti 35 + 777 2 T 2 / 35 - m 2 A 2j Bi 2 / + ^ 



(18) 
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It is also useful to examine the m — ► limit of the above expression as this case was calculated in Ref. 21] 
and can serve a valuable test for the above expression. By inspecting the occurring integrals we find the massless 
limit to correspond to 



ii i m __o 



ig 4 C A T F 
4(n - 1) 



n 7 



n— 
3 



+<T 4n - 6 



n 
n 



+ ^i) Tl35 + 16 if ~ q2 T) Tl235 + q2 (2n " 4) Bl2Bi5 



n 



1 



^2345 



(19) 



These terms are also, as expected, the only ones contributing to the gluon wave function renormalization 
constant. In other words, all divergent parts of the two and one loop integrals which vanish in the massless limit 
in the expression [l^ add up to zero identically. This in itself is an important check of the overall expression. In 
the heavy quark limit we can neglect the timelike component of the four momentum transfer q, i.e. qo = as was 
already mentioned before. This means that we do not need the longitudinal component of A4 gsei , however, we list 
it here for completeness: 
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ig A CATF 
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8 20 
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20 
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n- 
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A 2 B\ 2 i + 2^ 135 



2T< 
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28* 



^135 — 16T235 



n 



n^m 2 T 12 '3 5 - nT 235 + rJ-T nb + 4nA 2 B 45 - n^A 2 Bi 2 i + 8m 2 T 2U5 + im 2 T 2355 - 4Ti 35 



W (nT 23A5 + 2n^ 2 C 4 55 - 4m 2 T 23 455 - 4A 2 C 455 ) + ^ {4m 2 (T 23 5 - T 135 ) 
- (m 4 T 12 / 35 - m 2 Ti 35 + m 2 T 2 / 35 - m 2 A 2 B 12 > + A\ 



(20) 



A good check on the consistency of the employed decomposition is given by the absence of infra-red divergences. 
None of the two point amplitudes in this work is infra-red divergent to begin with, however, in intermediate steps 
of the calculation those do occur. An example is given above by the two integrals T 2 355 and T 23A ^ for which 
only the combination q 2 T 23 ^ — T 23 ^ is infra-red finite and this is how they enter into Eqs. 18 and 20. The 
function A 2 C^ only seems to have an infra-red divergence, however, in dimensional regularization it can actually 
be written as an ultra-violet divergence. This is done in appendix [B]. 

For the two diagrams that have an Abelian topology, Eqs. ^ and [?], we also give explicit results as usually 
only their sum is given in the literature [0, 10]. Here, however, we need both contributions separately due to 



the different color factors. In addition, Abelian and non-Abelian terms are separately gauge invariant and might 



display a different threshold behavior [14]. We find: 
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and for the longitudinal component: 
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Similarly, for Eq. ^ we get the following result: 
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(24) 



It can easily be seen that both parts of the two functions in Eqs. 21 and ^ multiplying are identical up 



to a minus sign when 23 is multiplied by the multiplicity factor 2. This is required by the gauge structure of the 
gluon propagator. Also their longitudinal parts add up to zero for the terms proportional to C F only. This just 
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checks the well known properties of the Abelian theory. It does not hold for the Ca parts of Eqs. 18 and ^ as 
they would get modified by the additional diagrams. These, however, were calculated in this work without the 
above reduction scheme as follows: 

We use the result of the integrated fermion loop which reads (omitting color and coupling constant factors) 



d n l Tr{- ffl {l/-fc + m)~f u {l/+m)} 
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(26) 



where rj is given in Eq. 85. For completeness, we also list the sum of the gluon and ghost contributions in the 



Feynman gauge [17], |36[ to the gluon propagator: 



7T, 
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87r 2r(4-e)(l-|) 
where £ is given in Eq. 85. Now W6 get the following result for J^A,gse^ 
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All the necessary integrals are given in appendix [B|. For the vertex correction graphs we arrive at the following 
representations: 
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where a is given by Eq. 126, p by Eq. 127 and a by Eq. 85. The remaining abbreviations read: 



B = -24(1 - x) + 8 + e(12(l - x) - 4) 
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(33) 
(34) 
(35) 

(36) 



The HQET Feynman rules of Fig. [I] project all three Lorentz indices to zero for A4 VCl . The completely 
antisymmetric nature of the three gluon vertex then implies that there is no divergence coming out of the internal 



fermion loop. Although Eq. 32 appears to possess a double pole, the " divergence" is actually finite when 
integrated over all Feynman parameters. We checked this directly with VEGAS and it indeed gives a well 
converged numerical answer. As for A4 VC2 , we integrate out the fermion loop as before, which yields: 
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-ig 6 C F C A T F 



VC-2 



2q 2 



d n k 
(4vr)" 



ir(k 2 ,m 2 



k 2 + 2kq 



(k + q) 2 k* 
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1 



(37) 



2(n-l) \(k + q) 2 k 2 2\{k + q) 2 k A q 2 k 2 q 2 (k + q) 2 J _ 
All the integrals left are given in appendix ^. 

3.1 Infra- Red Cancellations 

In this section we turn to diagrams which give integrals already present in an Abelian theory, however, which do 
not contribute in QED due to a cancellation that fails in the case of QCD. The reason is as follows: The color 
factors for the ladder diagrams are proportional to C F for the straight and Cp — Ca ^ f for the crossed ladder 



graph. The same structure is also present in graphs Eq. 10 and 12. In the sum of all four occurring ladder 
diagrams with one fermion loop plus M. VC3 and Moivc-, all terms proportional to Cp give a contribution that is 
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equal to the product of the one loop fermion graph with the Born contribution. This is an explicit example of 
the aforementioned exponentiation that occurs on the level of the potential. In an Abelian theory one thus has to 
omit these contributions. 

On the other hand, in QCD, we need to calculate the crossed ladder terms and keep only the — Ca £ f part of 
the above color factors. 

From direct inspection it is furthermore obvious that these diagrams contain infra-red (IR) divergent terms 
which have to cancel in the potential. It has been shown in Refs. [41, 21] that the sum of M c i + M V c 3 + M i vc 



is IR-finite. This requirement poses a strong check on the calculation and necessitates the calculation of the IR- 
divergent parts of a diagram that vanishes in dimensional regularization (A4 lvc), i- e - when UV- and IR-divergences 
are not separated. 

The presence of the square of the heavy quark propagator complicates the calculation of the crossed ladder 
diagram considerably as it makes the analytical separation of the double and single pole terms extremely difficult. 
We therefore found it most convenient to introduce a gluon mass A as an IR-regulator. This allows to explicitly 
differentiate between UV- and IR-divergences and provides a strong numerical check on the sum of all IR-divergent 
contributions. In this case we get the following integral representations for the unrenormalized and IR-regulated 
amplitudes: 



ig 6 C F C A T F f d n k 7r(k 2 ,m 2 )(k 2 -k 2 ) 

d 2 J (ATr) n (k + ie) 2 (k 2 -X 2 + ie) 2 ((k + q) 2 -X 2 + ie) [ ' 



-ig 6 C F C A T F f d n k 7r{k 2 , m 2 ){k 2 - jfeg) 
2a 2 J (4vr)« (k + ie) 2 (k 2 - A 2 + ie) 2 



-ig 6 C F C A T F . 2 o, f d n k 1 

M olvc = 7T-A TrO? > m ) / 7a wTTr 1 ■ no/? 2 \ 2 , ■ N ( 40 ) 

2g 4 J (471-)™ (ko + ie) z (k z - X z + ie) 

For the contributions of graphs ^ and 39 in which the kg terms in the numerator cancel the heavy quark 
propagator, no gluon mass regulator in needed. The sum of these feo-independent parts of M. c [ and M VC3 are 



separately IR-finite and indeed proportional to the integral 135 in appendix [B|. We therefore restrict ourselves 
to a discussion of the /co-dependent contributions only. In these integrals the i-e prescription is crucial in order 
to arrive at the correct location of poles and branch cuts the in the complex ko — plane. The presence of the 
fermion mass brings about a complicated integral over a general power, which in turn leads to a branch cut in 
the upper half of the plane. After integrating over ko in such a fashion one is left with an Euclidean integral over 



(n — l)-dimensions. More details and complete results are given in appendix B.l 
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4 Renormalization 



In Fig. ^| we list the relevant counterterms for the two loop diagrams of Fig. ||| and Fig. ||. The counterterms 
themselves contain non-local contributions, i.e. non-polynomial in the momentum transfer q, that have to cancel 
the non-local terms from the original amplitudes. The construction of the local wave function renormalization 
constants provides a powerful test of the correctness of the results presented both in section || and the appendices 
as they must combine successfully to arrive at the required local double and single pole terms. It might be 
helpful to expound on the general treatment of masses within the corresponding integrals and counterterms in 



the MS-renormalization scheme y, 22]. In the counterterm approach, their contribution is restricted to finite 
changes through the counterterms as the wave function renormalization constants are independent of the fermion 
masses. In other words, all pole terms that contain masses represent non-local infinities which must cancel in the 
sum of graphs contributing to the overall field strength renormalization. There is therefore no difference in the 
formal treatment of the mass parameter in graph ^ and any other graph. This is another way of saying that the 
parameters of a MS-renormalized theory are not physical. Rather, they are related to measurable quatities by a 
perturbative series in the physical parameters. 

We begin by presenting the results for the counterterms corresponding to Fig. ||. All two point counterterms 
correspond to the transverse parts of the gluon self energy contribution only, as these are the only relevant ones 
for this work. The graph M. gsei has two counterterms, one stemming from the fermion loop divergence (TI c t 1/ ) 
and one from the loop around the three gluon vertex (n c t 1 ). They are given in the MS-renormalization scheme: 



n 



X C A T F 



at i 



f 3(47r) 4 e J 



dx 



(8-6n)- -^F -l +2 



T x(l-x)\ -gVr 



T x(l — x) 



(An - 6) ( ^r(-l + f - x)) 1 - ?Vr(0(-^*(l - *)) 2 ) + (? 2 ("2 + 2n) • 



(0( A 1 



n — 1 



(41) 



Qig A C A T F f 1 



(4vr) 4 e Jo 



dx 



{An - 12) 



-m 2 r -1 + - 



T x(l - x) + 1 + q 2 x 2 T 



- (An- A)(m 2 + q 2 x)T( - --^x(l - x) + 1 - 4m 2 r -1 + - -^x(l - x) + 1 



2j\ m 2 I V 2/1 m 2 



T x(l - x)+ 
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vc 2 - ct: 



C7\ 



X 



f 




cl-ct: 




Figure 4: The two loop counterterms corresponding to the diagrams in Figs. [2] and ||[ Adding these contributions 

2 

to the original graphs removes all non-local functions from the occurring pole terms. The only exception are 
terms in the two point functions which only cancel in the sum of all two point diagrams as explained in the text. 
The fact that the tadpole diagram has no counterterm is already indicative of this cancellation. 
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gVr(-)l -±- x (l- x ) + l 



m z 



For the counterterm for M. gse2 we find: 



n — 1 



(42) 



±ig 4 C A T F q 2 



3n 



r 2 (1 



, 3\ r(i + |)r(i- §; 
+ in -2) T(l ' 



(43) 



ct2 3(4vr) 4 e(n - 1) \\2 "'"J T (2 - e) ' V" 2/ r (1 - e) 

where r\ and ( are defined in appendix [B| For M. gsei there is no counterterm as the subdivergence is independent 
of the mass which means that in dimensional regularization all the remaining integrals vanish. 
The pole terms for the respective terms, expanded up to O (e ) , thus read: 



ui + n 



O(e0) 
O{e ) 
O(e ) 



ig A C A T F q 2 
(4vr) 4 

ig 4 C A T F q 2 f 
ig 4 C A T F 18m 2 



_L i^. _ 
Qe 2 + 108^ ~ 

44 25 
~9e 2 + 27~e + 



3m 2 \ 

15m 2 
q 2 e 



(44) 
(45) 
(46) 



(4^) 4 e 

2 

These equations contain no non-local terms other than the — terms, which then have to vanish in the sum of 
all contributions to the non-Abelian part of the gluon wave function renormalization constant. Because of the very 
involved nature of the occurring non-local terms, this is already powerful evidence of the correct evaluation of both 
the two loop integrals as well as the decomposition of graph M gsei . Multiplying each graph with its respective 
multiplicity we find in the MS-scheme: 



{4[n t 1 + n ctl/ +n rfl9 



+ 2 



uf + n 



+ 



n 



}o(e0) 



ig*C A T F q 2 f 28 
V 3e 2 



71 

97 



(47) 



ct2 l ' r*J J ©(«<>) (4vr) 4 
This is completely local and thus demonstrates that the renormalization has been carried out properly and that 
the integrals given are correct. In order to further check this term we need the pole term from the "overlapping" 
Abelian two point diagram from Eq. |6| (which in QCD develops a color factor proportional to (C F — \C A ) in 
order to get the fermionic part of the overall gluon wave function renormalization constant Z3. The counterterm 
for M g seA reads 



IT, 



1 cr 



cti ~ e(4vr) 2 
and gives in agreement with [12]: 



{[ n ' + n ^Wo) 



C F -^f)T F q 2 ,{q 2 ,m 2 ) 



ig A (C F - %)T F q 2 ( 16 52 



(48) 



(47T) 



3e 2 9e 



(49) 
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Adding Eqs. [IT] and the Ca term of ^ gives the correct non-Abelian fermionic part of the gluon wave function 



renormalization constant ((times Jj) see Ref. [36] for example) in the Feynman gauge: 



Z c - = 9 " CaTf - l\ (50) 

3 fermionic (471")^ \ 3g^ 6/ 

This testifies to the overall correctness of both the decompositions used as well as all the integrals listed in the 
appendices! 

For completeness we also give the counterterm for M gse5 , which in the MS-scheme must be treated in the same 
way as the graphs before. All divergent terms proportional to m 2 cancel the corresponding non-local infinities in 
Eq. M. 



n cts 



with 



4 i g 4 C F T F 
e(n- l)(4vr) 2 

-2Am 2 B 22 - (V + 16m 2 )£i2 + 8^2 



n ( \2m l B 22 + 2q 2 B 12 - AA 2 - 12q z m 2 C 122 ) + (2Aq 2 m 2 - 48m 4 )Ci 2 2 



(51) 



( [n? + n cte j } 

U * 5 JIo(eO) 



ig A C F T F q 2 
(4vr) 4 



8 8 
+ 9e 



(52) 



It is an important difference to the massless case that the counterterms |48| (rather its C F part) and 51 are 
not related by a simple minus sign as implied by the Ward identity. There is an additional constant term 4m 
which gives new contributions. For the purely Abelian fermionic part of the gluon wave function renormalization 



constant in the Feynman gauge we find in agreement with Ref. [36|: 



z: 



C F 



g C F T F (_4 



(53) 



^fermionic ^471"^^ 

The cancellation of the higher order (double) pole is a characteristic feature in QED that holds to all orders 



For A4 VCl we do not need to remove non-local terms as the fermion loop is finite due to the projection of all 
three Lorentz indices to zero. It is easy to check this by calculating all divergent pieces after the integration of 
the fermion loop. All that is left is the divergence from the remaining integral which has to be subtracted in the 
usual MS-fashion. This is indicated in Fig. 0. The explicit pole term is given by: 



[M 



vcl 



ig 6 C F C A T F 



0(e°) 



(47T) 



4„2 



(54) 
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in agreement with the massless case 2C]. In the case of M VC2 we do have non-local terms, and the counterterm 
reads: 

Aig 6 C F C A T F 7?! /-i / nvT (§) , (1 + v)T (l + f ) , 



M*2.* = n s / -o — r + 



3(47r)Ve Vo^V 2(^«(l-«))f ' (z£t,(l-«))3 



r(-) _ | n,r(f) _ (l + tQrft+j) \ 



2(n - 1)(^(1 - 8(n - - 4(n - - «))* 

Adding Eq. ^ with the appropriate normalization and color factors to the result given in |37] does indeed give 
completely local double and single pole terms as required in dimensional regularization after the subdivergences 
are subtracted: 

ika a d i ig G C F C A T F ( 1 5 \ 

It demonstrates that indeed all non-local divergences are canceled and agrees furthermore with the pole terms 



obtained in the massless analysis [20]! It should be noted that all the integrals needed were already used in the 
Aigse-2 calculation for which such a strong internal consistency check was performed just above. All the required 
expansions above were carried out with the help of MAPLE in face of the complexity involved. As mentioned 
before, also the translation into FORTRAN was handled by MAPLE as to reduce possible accidental errors. . 

4.1 Counterterms with Gluon Mass 

At this point we need the counterterms of the IR-divergent contributions, Aid, M VC3 and M. \vc- As indicated 



above and expressed in Eqs. 38, ^ and 40, these were regulated by introducing a gluon mass regulator. The 
remaining UV-divergences are treated as above in the context of dimensional regularization. We therefore have 
to calculate all counterterm contributions that occur for gluon propagators with a gluon mass. Without such a 
dimensionful quantity, only the crossed ladder diagram would yield a counterterm in dimensional regularization. 



We again use the gluon mass only for /^-dependent terms as explained in section 3.1. This is indicated below 



The results are obtained in a similar way as for the corresponding amplitudes, first integrating over the heavy 
quark propagator in the complex fco-plane with a subsequential (n — l)-dimensional Euclidean integral remaining. 
The results are obtained straightforwardly as there are only pole terms and no branch cuts in the counterterm 
contributions. We find for the gluon mass regulated terms: 
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Figure 5: The sum of the A 2 -dependent amplitudes and counterterms M k c ° + M k v % z + M ivc + + ^v°c 3ct - 

Circles correspond to a choice of q 2 = —WGeV 2 and m = m c , triangles to q 2 = — WOGeV 2 and m = m c while the 
lower curve (squares) has q 2 = — WOGeV 2 and m = mi,. The overall normalization neglects color factors and the 
coupling strength. All data are obtained by using 10 6 evaluations per iteration with VEGAS and 100 iterations. 
The statistical error is indicated and smaller than the symbols where invisible. The sum for each of the displayed 
sets of parameters is clearly independent of the IR-gluon mass regulator A as expected. 



LI + + ♦ * • 



"A A A A A A A A 



, _ tfWJU-rluj) f' dv _ 2 ■ (57) 
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Mk0 = gjg^uWr^ 



3^ 2 (4vr)ir(|) (£) 

For completeness we also list the remaining counterterm stemming from the fco-independent part of M c f. 



(58) 



= m^4Mr (lt |)r[ f ) 
*' 3^(4,)^r(i- ! )( 1 #) 1+ - 

The gluon mass terms that occur in the expansion of the original as well as the counterterms above in powers 
of e in the pole terms of dimensional regularization represent now non-local divergences which have to cancel 
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in the same way as terms containing m 2 or non-polynomial functions of q 2 . The remaining IR-divergent pole 
terms are contained in the form of logarithmic divergences in A. Fig. || demonstrates that in the sum of the IR- 
divergent amplitudes plus their corresponding counterterms no A-dependence is left within the statistical errors. 
For convenience, three sets of values for q 2 and m 2 are displayed while the renormalization scale remains fixed. 
We have checked that it also holds for a variety of other choices of parameters. Some need fewer evaluations to 
converge while others need up to 10 7 per iteration. 

It is perhaps interesting to note that the crossed ladder diagram, naively only singly IR-divergent, actually 
possesses a quadratic divergence in log(X) which cancels the (also unexpected) quadratic divergence in the Abelian 
vertex correction term. The remaining UV-divergent pole terms in the MS-scheme are found to be: 



M k J+M k c ? 



M ko +M' 

J l vc 3 ~ J l vc 3ct 



ffco 



+ M k clc 



O(e ) 
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(4tt)V ( 
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ig 6 C F C A T F . 


f 16 


16 

+ 97 


0(e°) 


(4tt)V ( 


v"3? 



(60) 
(61) 
(62) 



which states that the counterterm in case of MJ completely remove all pole terms in e. It is also clear that 
all non-local terms are removed by the appropriate counterterms as was expected. In order to compare this with 
the results obtained in the massless case one would need to differentiate between ejjv and em- 



5 Numerical Results 



At this point we have calculated all diagrams that contribute to the massive fermionic corrections to the heavy 
quark potential that were previously unknown. In the previous section we demonstrated that the counterterms 
successfully remove all non-local divergences and that the MS-subtr action terms coincide with the massless limit. 
The complexity of the explicit results given in the appendices raises some questions about how stable a numerical 
integration over up to four Feynman parameters is with VEGAS as well as about the correctness of the finite 
terms of these expressions. An ideal test is provided by the results obtained in Ref. [^] for the massless limit. 



Fig. H contains the results of the IR-finite two loop amplitudes from Figs. || and |3| in section |2[ The tadpole 
diagram vanishes trivially in that limit so that only the six graphs shown remain. The sames choices for q 2 
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Figure 6: A comparison of the six amplitudes M. gsei , M gse2 , M. gse ^ M. gS e 5 , -M. VC1 and M VC2 with the massless 
limit (dashed lines) |2(| in the MS-scheme. Solid circles correspond to a choice of q 2 = —1.5 GeV 2 , open ones to 
q 2 = —4.5 GeV 2 . fi = 0.31 GeV in each case. Each graph begins to deviate from the massless limit only when m 2 
is of the same order as — q 2 as expected. These results were obtained after 10 6 evaluations per iteration and after 
50 iterations. The statistical error is smaller than the size of the symbols and the normalization neglects color 
factors and the coupling strength. 
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Figure 7: A comparison of the sum of amplitudes M c i+A4 VC3 +M ivc plus their MS-counterterms with the massless 
limit (dashed lines) pCfl . Solid circles correspond to a choice of q 2 = —1.5 GeV 2 , open ones to q 2 = —4.5 GeV 2 . 
fi = 0.031 GeV and A 2 = 10 -8 in each case. The sum begins to deviate from the massless limit only when m 2 is of 
the same order as — q 2 as was the case for the other graphs. These results were obtained after 10 6 evaluations per 
iteration and after 100 iterations. The statistical error is smaller than the size of the symbols and the normalization 
neglects color factors and the coupling strength. 



and the renormalization scale /i were made in all six plots. Since the results of Ref. pi were calculated in the 



MS'-renormalization scheme, we use 

»MS = \[^^ms ( 63 ) 

It is clear from these results that deviations from the massless limit only occur when m 2 ~ — q 2 or greater. 
This was of course expected and the motivation for this calculation. A similar dependence is observed for the sum 
of the three IR-divergent amplitudes from Fig. ^ in section ^. Here it is impossible to compare on an amplitude 
by amplitude level since a different IR-regulator was used. Only the sum of infra-red finite contributions can be 
compared at the two loop level. We checked explicitly that by replacing log(X) with |, only the double pole terms 
can be seen to be identical. 
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The single pole terms differ and so do the finite contributions for each amplitude. In the sum, however, the 
IR-divergent pieces cancel (as demonstrated in Fig. ||), and here we can find a meaningful comparison. Fig. [?] 
demonstrates that the correct massless limit is indeed recovered. The numerical accuracy in terms of the statistical 
error from the VEGAS Monte Carlo integration is actually included in the figures. It is better than 1% though, 
and thus not visible on the scale of the plots. The gluon mass regulated graphs were evaluated over twice as many 
iterations (100) compared to the graphs from Fig. |6| as the required cancellations are numerically more unstable. 
In both cases the number of evaluations per iteration is 10 6 . 

Fig. H displays the sums of all non-Abelian as well as the sum of all Abelian fermionic contributions to the 
heavy quark potential. In addition we included the one loop corrections (bottom) in the MS-scheme (omitting 
coupling constants) as given in Eq. |26|. The simple logarithmic behavior of the massless one loop result is clearly 
visible and asymptotically approached by the massive curves. The sign of the one loop correction is opposite to 
the two loop Abelian result, reflecting the fact that effectively for large momenta (3q ED — > (p*Q ED ^j (in the 
massless case, with /3q ED = — |). The relative size of the mass effects are comparable for the one and two loop 
corrections. 

The massless two loop results can be seen to possess the expected double logarithmic contributions. The 

2 

massive two loop results show an almost completely opposite behavior for low values of ^-j. At the flavor 
thresholds, though, both contributions increase the value obtained from the massless case by the same (relative) 
order of magnitude. The overall corrections are much larger in absolute terms for the non-Abelian case, partially 
due to an extra factor of Ca, while in relative terms the Abelian corrections are bigger. In the high energy regime 
both graphs show that the massless limit is approached asymptotically. 

The complete massive fermionic two loop contributions to the heavy quark potential are presented in Fig. |9[ It 
can be seen that the overall curve is dominated by the non-Abelian threshold behavior (partially due to the extra 
factor of Ca)- The "m c -graph" (triangles) matches the massless case for lower values of —q 2 as m 2 <C m 2 . At the 
respective thresholds we find roughly a 33 % deviation relative to the massless case. This could be very significant 
for applications where quark masses are expected to play an important part. Furthermore, the physically defined 
effective charge ay(q 2 ,rn 2 ) incorporates quark masses naturally at the flavor thresholds and is also analytic. 
Thus, there is no problem of evolving the strong coupling constant through these thresholds and one never needs 
to impose matching conditions. At high values of q 2 the theory becomes massless and reproduces the leading 
logarithmic terms obtained by the /3-function analysis as these coefficients are scheme independent through two 
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Figure 8: The mass dependence of ay at one (bottom) and two loops. The two loop case is displayed in terms of all 
Abelian terms (left) and for all non-Abelian terms (proportional to Ca)- Triangles denote m 2 = m 2 = (1.5GeV) 2 
and open circles m 2 = m 2 = (4.5GeV) 2 . The massless case is also included (line). All curves have the same value 
of the renormalization scale jJL = 0.031. It is clearly visible that the flavor threshold behavior is quite similar in the 
three figures with an opposite tendency for low values of — q 2 in the two loop case though. The one loop corrections 
have an opposite sign relative to the Abelian two loop corrections. The coupling constants are omitted. All cases 



approach the massless limit when <C 1. 
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Figure 9: The complete two loop mass dependence of dy = — — for m 2 = m 2 = (1.5GeV) 2 (triangles) and 
m 2 = m 2 = (4:.5GeV) 2 (open circles). The massless case is also included (line). In all three curves we use 
[i = 0.031. The deviation from the massless case at the flavor thresholds is of order of 33% and is dominated by 
the new non-Abelian contributions 

loops in a massless theory. 

The above analysis can also be helpful for the incorporation of massive fermions in lattice analyses as the 
heavy quark potential is defined by the gauge invariant vacuum expectation value of the Wilson loop in Eq. [|. 
For a direct application of the presented results, a recently proposed way of incorporating quark flavor thresholds 
by relating the "natural" heavy quark potential m^-dependence to an effective continuous and smooth function 
nF(—q 2 ,m 2 ) |k| seems to be a promising candidate. 

6 Conclusions 

We have calculated all the necessary integrals for the non-Abelian massive fermionic corrections to the heavy quark 
potential through two loops. They describe the analytic flavor thresholds of the physical coupling ay(q 2 ,m 2 ). 
The new results were obtained by using a mixed analytical, computer-algebraic as well as numerical approach 
and strong consistency checks were performed by observing that all non-local divergences cancel by adding the 
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appropriate counterterms. In case of the complicated two point diagrams it is found that the weighted sum of 
all diagrams gives the correct local gluon wave function renormalization constant. The renormalization constants 
were given explicitly. 

It was also checked that no spurious infra-red divergences were introduced by the implemented reduction 
scheme as they are present in the intermediate steps of the calculation. For the explicitly IR-divergent diagrams 
we proved that no physical results depend on the introduction of the gluon mass regulator A. This is a consequence 
of the color singlet state of the external heavy quark sources. 

All physically interesting and gauge invariant finite parts were integrated with VEGAS pf| and found to agree 
perfectly with the massless results of Ref. [^2| in that limit which actually checks this part of the analysis in |2l[ . 
The difference to the massless case around the charm and bottom flavor thresholds was found to be roughly 33%. 
The size of this effect can have important consequences for processes in which one cannot neglect these masses as 
well as for the evolution of the strong coupling constant through analytic flavor thresholds. 



A Decomposition of Two Loop Tensor Integrals 

For the gluon self energy graphs A4 gsei , A4 gse4 and Ai gS e 5 we chose to not do the fermion loop integral first, 
as we did for all vacuum polarization insertions, but to decompose the occurring tensor integrals into a linear 



combination of scalar two loop integrals. The scalar integrals entering in the expression given in Eq. [1^ (or 2C) 
will be treated in detail in the next section together with all other integrals needed in this work. 

We work in n space-time dimensions, n = 4 — e, and for the two loop integrals we use the following notation: 



[1] = {l + qf - m 2 , [2] = l 2 - m 2 , [3] = (I - kf - m 2 , [4] = (k + q) 2 , [5] = k 2 (64) 

I denotes the loop momentum of the massive fermion loop and k the remaining loop momentum. A prime like 
[2'] = I 2 denotes the massless propagator with the same momentum as the unprimed. While there are different 



possible technical approaches to our desired decomposition, such as the one recently suggested in Ref. [24], the 
general method we use follows that of Ref. . We also denote integrals with squares of "denominator" terms in 
the numerator u Y n - integrals and pure two loop scalar integrals by "T". 

In the following we use various symmetries between Y- as well as between T-integrals. For instance 



Y, 



1345 



Y 1 
*2345 



^134 



'-235 



(65) 



For two loop scalar integrals that are actually a product of scalar one loop integrals we use the respective one 



loop notation of Ref. [|ll|] . All of the decompositions were programmed in FORM [13] and lead to the following 
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7n 



+m Z T 235 + -g 4 r 2 345 + 77^26455 ~ m 2 g T 23 455 ~ q^A 2 C^ + 7-?™V?12'35 + ■-77VT135 



n 



2„2r 



6 



12 



77 



<? Z 7i35 + nq z A 2 B i5 - '-q z 'A 2 B 12 , + 2m 2 q 2 T 2U5 + m z q z T 2355 - q z T 135 
b 



,2„2r 



(66) 



2345 



1235 



2 1245 



231 



Y 



135 



r 235 



245 



(2vr)" 



^2^45 + 



d n k 



d n l 



li 2e {{l + qf -m 2 ) 



(2tt)' 



(27r)« (P-m 2 )((/-A:) 2 

^235 — ^135 + 9 2 ^2345 



m 2 ){k + q) 2 k 2 
fi 2e (k + q) 2 



(2Tr) n ({I + q) 2 - m 2 ){l 2 - m 2 )((l - k) 2 - m 2 )k 2 
A 2 B 12 + g 2 Ti235 + 7j ( T 2'35 + -42 £22' - ^2^12' - T 135 - q 2 A 2 C 122 > + {m 2 - q 2 )T 12 > 

^((l - k) 2 



d n k 



d n l 



m 



(2vr)™ J (2vr)" ((I + q) 2 - m 2 )(l 2 - m 2 )(k + q) 2 k 2 



A 2 B, 



2-D45 



--B12-B45 



/^ 2e ((Z + g) 2 -m 2 ) 



d"Jfe 

(2^7 (27r)"(/ 2 -m 2 )((/-A:) 2 -m 2 )(A; + g) 2 



^ 2 + m 2 



d n k 



T 2 '35 

d n l 



A 2 B\ 2 i — T135 + m Ti2'35 
^{k + qf 



+ q 2 (T 135 +A 2 B 12 , -m 2 T] 



12'35 



(2vr) n 7 (2vr)™((Z + g) 2 -m 2 )((/-A;) ; 



777/ 



! )fc 2 



~ (A 2 + g 2 T 13 5 



// 2e ((/ + ^) 2 -m 2 ) 



l 2 

d n k 
(2vr)™ 

(2vr)" 7 (2vr)" (Z 2 - m 2 ){k + q) 2 k 2 



m 2 (T 2 > 35 - A 2 B 12 i - T135 + m 2 r 12 '35 ) + q 2 [A 2 B 12 > - m 2 T : 
2 



12'35 



(2vr)" (I 2 - m 2 ){(l - k) 2 - m 2 )k 2 



q 2 T 235 



d n l H 2e {(l + qY 



m 



q 2 A 2 B i5 



(67) 

(68) 
(69) 

(70) 

(71) 
(72) 

(73) 



For the remaining two diagrams, A4 gsei and A4 gses , we have slightly different denominators as is evident from 
Eqs. ^ and |?]. It is possible, though, to relate these to the conventions given in ^34] with the exception of the 
finite scalar integral T^ 345 which is given in Eq. |128| . "^4" denotes the fact that the topology of these diagrams is 
Abelian. Below we list the ^-functions we need for the required decomposition with terms on the l.h.s. having the 
denominators of the original integrals and given in terms of functions on the r.h.s. which are using the conventions 
of Eq. |4: 
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A v l 
-^2345 


— A 2 Bi2 + T 2 35 — T135 + g 2 Ti235 — y"!235 


(74) 


A v 3 
1 1245 


= 2A 2 B 12 + (2m 2 - \)B 2 2 


(75) 


A V A 

ZOOO 


= ?235 + 9 2 ?2235 


(76) 


A V 2 
1 135 


_ iy5 _ iv4 _ Avl _ v 1 
— 2 234 — M35 ~~ 2 234 _ 2 234 


(77) 


2 235 


= 9 2 ?235 


(78) 


1 235 


= Q 2 T 2 35 


(79) 


2 255 


= A\ + q 2 A 2 B 22 


(80) 



B Two-Loop Integrals 

In this appendix we give the explicit results for all the integrals needed in the calculation of the two loop fermionic 
corrections to the heavy quark potential. These include all the scalar two loop integrals occurring in the decom- 
position of the gluon self energy graph Ai gsei in section ^ as well as the remaining tensor integrals needed for the 
remaining contributions. Since the potential between two infinitely heavy color test charges represents a physical 
quantity, all integrals presented are real due to the spacelike value of the physical momentum transfer q 2 < 0. For 
this reason we found it convenient to adopt both analytical as well as numerical methods for the implementation 
into FORTRAN. Wherever possible we proceed with the integration of the remaining Feynman parameter integrals 



and where this becomes too involved, we integrate the remainder with the Monte Carlo integrator VEGAS, [25]. 
The notation is as follows: 

The following Feynman parameter identities 1 26 ] are very useful and were employed in all integrals in this work: 
r(m)L...k-n : 1 \ • m ~ 2 , : jz (81) 



a\...a m Jo Jo (aiui...u m -i + a 2 ui...u m - 2 {l - u m -i) + ... + a m (l - ui)) r 

1 _ T(a + m rt -t' (82) 



a a b? T(a)T((3) J (au + 6(1 - u)) a+ P 

i _ r (a+/3+7) 4,4 (■r;y-.)r(i- ; .r; m 



a a bPci Y{a)T{p)T(j) Jo Jo {auv + bu(l - v) + c(l - ii)) a+/3+7 

1 r(a + /3 + 7 + <5) f\ 2 f\ f\ {uvw) a - 1 {uv{l-w)f- 1 {u{l-v)y- 1 {l-u) s - 1 foA . 



a al,p &d s r(a)r(/3)r(7)r((5)7o Jo Jo (auvw + buv(l - w) + cu(l - v) + d(l - u)) a +' 3 +~f+ 5 
We use the following abbreviations in addition: 



Airfi 2 



a 



u + (1 — u)x(l 



a 



u + (l-u)(l 



(85) 
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A s 4W';»)'_!!M) t ^ (86) 

m z \ or a I a 

~ _ q 2 (u z {l-v) z u(l-v) \ I 



-,2 



a 2 


a 


u 2 {l-v) 2 


u(l — v) 


a 2 


a 


u 2 (l-v) 2 


u(l — v) 


a 2 


a 



+ = (87) 



A' = -^1 ■ ^^1 + 



\ or a a 



where n is the dimensional-regularization mass parameter [23]. All results are given in terms of their dependence 
on e and would have to be expanded with the factors given in the explicit results of section | up to O (e). The 
results in this paper were obtained by employing MAPLE to do the required expansion and are too cumbersome 
for explicit presentation. 

We start with results of the following simple scalar one and two loop functions: 



= <Tl = imW-l + f) 

1 ~ 1 (27r) n {I 2 - m 2 ) 16vr 2 K ' 

d n l // _ir ? §r(f) 

22 = ' (2vr)« {I 2 -m 2 ) 2 ~ 16tt 2 (90) 

2 - ~ 1 (2ir) n (I 2 - m 2 )l 2 16vr 2 (l-f) ( ' 

Bl2> = / 7 ^771 77T oTTo = / ^ i ^ — 92 

' (2^ {{I + q) 2 - m 2 )l 2 J ( 4vr )2(^ x (l_ x ) +x )f V ; 

= ' (27r)™((Z + (? ) 2 -m 2 )(Z 2 -m 2 ) " /„ X (4vr) 2 (^x(l - x) + 1)1 (93) 

^ ^ _ yi ^fr(f) _ iC fr(f)r 2 (i-f) 

" ' (27r)»(* + ?) 2 fc 2 io (4^) 2 (^(l-x))§ (4vr) 2 r (2 - e) 

= <ffc // < fr(-f)r(i-f)r(i + f) 

'" ' (27r)» (fc + q) 2 k 4 g 2 (47r) 2 r(l-e) 1 ' 

( 122 = ' (2^Y{(l + q) z -m 2 ){l 2 -m 2 ) 2 ~ ~ Jo * (4vr) 2 m 2 (^x(l - x) + l) 1+ § ^ 

, d n l ^ _yi r 1 ir ? fr(i + f)x-f 

l " = ' (2vr)«(a + ^) 2 -m 2 )(/ 2 -m 2 )/ 2 = io X io 2/ (4vr) 2 m 2 (4(x(l- 2 /) 2 -l + y) + l) 1+ f 1 j 



A very useful integral for 97 is given by 



1 , f 1 1 2 , -i / / « \ a ~ 1, / \ / \ 

ax / dy-——- -75 r — = — V a z — 4a tanh w /oo(l — a) (98) 

o Jo y (a(x(l-y) 2 -l + y) + l) a VV a - 4/ a yv ; y J 

This integral is needed in order to analytically separate the divergent pieces since C122' is multiplied by A 2 in 

the solution for Eq. |(| 
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_d»k_ r d n i ^ mVr(f)r(-i + 6)r 2 (i-§) 

~ i ' ~ ' (27r)"7 (2ir) n l 2 ({l-k) 2 -m 2 )k 2 {Att^T (2 - § ) 1 J 

d n fc /• d n / ^ e _mVr 2 (!)r(-i + 6)r(i-0 



[ 235 



(2vr)«7 (2tt)™ (Z 2 - m 2 )((/-£;) 2 -m 2 ) A; 2 (4vr) 4 r (2 - §) T (e) 

The reason why the following integrals cannot be given in such a simple form is the presence of the external 
momentum transfer q in addition to the masses. In order to extract the infinite pieces from the next integral T135, 
we repeatedly use the following propagator identity: 



1 1 2lq + q 2 



(I + q) 2 — m 2 I 2 — m 2 (I 2 — m 2 )((l + q) 2 — m 2 ) 

It then follows that 



(101) 



with 



d n k f d n l \i 2e 
Tl35= ' (2^7 (2vr)« ((/ + q) 2 - m 2 )((l - k) 2 - m 2 )k 2 - T ™+T a +T b + T c , (102) 



_ d n k f dH ^(2lq + q 2 ) g 2 T?er( e )r(e)r(l _. )r(l+ e ) 

"' ~ ' (2n) n J {2ir) n {I 2 -m 2 ) 2 {{l-k) 2 -m 2 )k 2 " (4^) 4 r (2 - § ) T (1 + e) 1 ' 



n 



d n k r d n i ^{2iq + q 2 ) 2 4gVr(f)r(e)r(i-f)r(i + §) 



(2tt)«7 (2Tr) n (l 2 -m 2 ) 3 ((l-k) 2 -m 2 )k 2 n(4vr) 4 r (2 - § ) T (1 + e) 

. g 2 (^ + ^)^(f)r(i + e)r(i-f)r(2 + f) 

2(4vr) 4 r (2 - f ) r (2 + e) ( j 

In passing we note that 

T 223 5 = ~\T a (105) 

The last term T c has only a simple pole in e which is, however, buried in the Feynman parameter integration. 
This is a quite common problem that arises because of the T-factors in Eqs. |82] and 83. We take "u" to be that 
Feynman parameter and for our purposes it suffices to write the following identity: 

1 du(l - u)^ l f(u) = -/(l) + f 1 du(l - u)!- 1 (/(«) - /(l)) (106) 

e Jo 

The respective terms for T c = - J ^ J ( g„ (fa _ ma)a((f ^^^j f 3 _ fc ^_ m ^ )fc ^ are: 
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and thus 



2 6 _ 12^1^1 

a -r(i + e) 



A 1+e 



(107) 



/(l) = f dx f dv 
Jo Jo 



2(4tt) 4 x2 
q 2 6-12(1- 




(108) 



Although this result for T135 is correct, it is numerically unstable in the massless limit because of terms of 
order which have to cancel as m 2 — ► 0. A way out of this calamity as well as a very good check on the 
correctness of our result for this integral is to use the propagator identity 101 for n^fqp instead after shifting the 
loop momenta. This yields 

T IM =T 2M -rt M5 -<< M || M >> (109) 

The result for T2345 is given below and the last term in the equation can easily be found to be 2q 2 u<yl ~ v ^ 
times the expressions for the scalar integral. This term just stems from the momentum shift k — ► k! — q u ^~ v ^ ■ 
Numerically, away from the singularity at m = 0, both solutions agree. 

In similar ways we treat the following more complicated integrals, always calling "it" the Feynman parameter 
that contains an additional divergence if /(it)-terms are quoted. The desired value for the respective integrals are 
understood to follow from an expansion in e of Eq. |106|. For 



- 2345 



d n k f d n l fi 



2c 



(2Tr) n J (2n) n (I 2 - m 2 )((l - k) 2 - m 2 )(k + q) 2 k 2 



(110) 



we get 



/(«) - - t dx f 1 dv ggg? , f(l) ^-f 1 dv ^ (111) 

Jo Jo (47r) 4 a 2+ 2A e JKJ Jo (4tt) 4 (-<v(1 - v)Y 



Similarly, 
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T _ f d n k r <ri £ 

1235 ~ J (27r)™7 (2tt)" ((/ + g) 2 - m 2 )(/ 2 - m 2 ){{l - k) 2 - m 2 )k 2 K ' 

with 

/(«) = - r <** f ,f r Ai^T £ , /(i) = - r <** r *> - u — — - ens) 

Jo Jo (4ir) 4 x2a 2+ 2A € Jo Jo (4tt) 4 x2 (-^v(1 - v) + l) e 

For 



T _ f d n k f <Pl ^ 

12 ' 35 7 (27rW (27r)™((/ + g) 2 -m 2 )/ 2 ((/-fe) 2 -m 2 )fc 2 1 j 



(2tt)" y (2tt)" ((/ + g) 2 - m 2 )/ 2 ((Z - fc) 2 - m 2 )fc 2 

we find 



f(u) = - [ dx f 1 dv r ^ V y )U € ^ e , /(l) = - f 1 dx f 1 dv 5 (115) 

Jo Jo (4vr) 4 x25 2 +2A' Jo Jo (4vr) 4 x2 (-|U(1 - v) + 1 - t>) e 

The infra-red finite integral 

/■ d n k f d n l ^{k 2 + 2kq) 
2455 ~ y (2vr)«y (2tt)» (/ 2 -m 2 )(k + q) 2 k* { ' 

is a product of two one loop functions which are given by 

. _ imVr(-i + f) 



(4vr) 2 

y (4tt) 2 \2(u(l-u))i \2J (n(l-u)) 1+ f V 2) J y ' 

and in dimensional regularization we have ^2455 = ^2-^455 = (f^C^, where ASzC denote the respective one 
loop scalar integrals. For the infra-red finite combination 



(117) 



T -„ 2 T T - f d ' lk [ d ' li ^(k 2 + 2kq) 

123455 = q T 234 55 - T 235 5 ~ ~ J 7^ J rp _ m 2 )((Z _ k) 2 _ m 2 )(k + q W (U9) 



d n k r d n l fi 2e (k z + 2kg) 

(2^ J (2vr)™ (I 2 - m 2 )((l - k) 2 - m 2 )(k + q) 2 k 4 

we get two u f(u)" terms, distinguished below by capital (containing double pole terms) and lower case (with 

only simple poles) letters: 



F(u) ee ['dxf'dv nriT } e jf. V , F(l) = f'dv nr]eT 2 ie)V (120) 

Jo Jo 2(47t) 4 q j+ 2A £ Jo 2(4tt) 4 (-^(1 - v)Y 
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f\ [ l j a r] e T(l + e)u 2 v f 1 (1 + «)t/T(1 + e) , 101 , 
/(«) = dx dvA ' s ; o.a A11 , , /(I) = - / ^ 7 A -r- ' 2 77" ( 121 ) 



JO 



"(47r) 4 a 3+ tA 1 + £ ' Vo (4tt) 4 (— - «)) 



., . it 2 (l — w) 2 u(l — v) . 



For 



12235 



2c 



(2vr)« 7 (2tt)" ((/ + q) 2 - m 2 ){l 2 - m 2 ) 2 {{l - k) 2 - m 2 )k 2 



we find 



/(«) 



lo (4vr) 4 x2a 3+ 2A 1 + e Jo Jo (47r) 4 xi(-Xr?;(l - v) + l) 1 

d n k r d n l n 2t 

(2n) n J (27r)« {i 2 -m 2 )((l+q) 2 ~m 2 )((l-k) 2 -m 2 )(k+q) 2 k 2 



The completely finite integral T 12U5 = I Mr f ( 2 tN^- m ^(i + rf-m4i-tf- m ^(H^^ is g iven b y ; 



(123) 



f 1 dx r 1 dv ^ T }\ + ;}fy , /(i) = C d X f 1 dv . ^ 1+e > (124) 

Jo Jo (4vr) 4 x2a 3+ 2A 1 + e Jo Jo (47r) 4 x5(-^(l - v) + l) 1+e 



m [ l i [ l i [ l i [ l i ?fr(l + e) xu(l-ii)i 

Ti 23 45 = / dx dy du dv ' v J — — i- — ^ 125 

Jo Jo Jo Jo m 2 (4vr) 4 a 3+f / <? 2 ^ 2 _ A , « 1+( 

e can of course be set to zero in the above expression and we use the following abbreviations: 



cr = u(l-u) + (l-ti)(l-i/)aj(l-aj) (126) 
p ee u(l-<t;) + (l-ii)(x(l-y)-x 2 (l-y) 2 ) (127) 

For the "Abelian" gluon self energy graph M qseA we need another completely finite integral with five denomi- 



nators, namely T& U5 = J^tJ ( ffi„ ( ;2_ m . )(( ; +g) ^ m 2 )(/ ^) 2 ((fc +9 )^ m 2 )(fc^m 2 ) • Here we find 



i /"I f 1 7?T(l + e) xu(l-u)i 



o Jo Jo Jo 



T\1 3 45 = /_ /. % /_ du I dv '^; A l^' 77' ,7 " g(1 _ u) 1+e ( 1^) 



m 2 (4vr) 4 a 3+| /V (g»_£\ + 

\ m 2 \a 2 a J 



Again, we can savely set e to zero like above. The following integrals are needed for the diagrams where we 
integrated out the fermion loop first, with -ir(k 2 ,m 2 ) taken from Eq. 



d n k fi e ir(k 2 ,m 2 ) f 1 f l m 2 Y (-1 + e) x(l - x)(l - n)"f rf 



(2tt)" k 2 J J 327r 4 a 2 -f 



1-e 



d n k n € 7r(k 2 m 2 ) f 1 f 1 m 2 r (-1 + e) x(l - x)(l - u)~§ rf (-^ ux ^ g) + A 

, f } = dx du — L 1 \ g! g Z (130) 

(2tt)" (& + g) 2 Jo Jo 32vr 4 a 2 -f V 7 

d-fc 2^vr(A ; 2 ,m 2 ) f\ f 1 -g 2 m 2 r (-1 + e) x{\ - x)u(l - n)"f rf (-^*&=& + l) ^ 

' ax / an — j (131) 



(2tt)" (k + q) 2 Jo Jo 167r 4 a 3 -2 
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Below we split again into f(u) and /(l) terms. For 



d n k /j, e ir(k 2 , m 2 ) 
(2vr)™ (k + q) 2 k 2 



(132) 



we find: 



1 rl T (e) ux(l — x)r] e 
o ""Jo "' 32vr 4 a 2+ §A e 



f(u) = - I dx / dv V J - v 2+ ? - f ' (133) 



For 



W io i 32vr 4 f~' ?M ^ ) 



d n k (k 2 + 2kq)n e ir(k 2 ,m 2 ) . . 

— (135) 



(2vr)™ (jfe + g) 2 & 4 

there are two contributions corresponding to terms with double poles (F) and only single poles (/): 



. f 1 , f 1 , r (e) nu 2 vx(l — x)n e 

F(u) = - dx dv—-t ^ J — 136) 

V ' Jo Jo 64vr 4 a 3+ 2A e V ; 

F(l) = - ['dx / V (e) 7*i\- g ff (137) 

,l rl q 2 T (1 + e) n 2 ra(l - x)rf ( - 2^±) 

f(u) = dx dv-- ° ^-L (138) 

' Jo Jo 32m 2 7r 4 a 3+ fA 1 + e V ' 

Jo Jo 32?r 4 ( q 2l ) 

B.l Two Loop Integrals with Gluon Mass 

In this appendix we give details about the evaluation of the IR-divergent integrals of section 3.1. The contributions 



containing heavy quark propagator terms were regulated using a gluon mass regulator and lead to the following 
general integral over fco : 

r = r ^° I rum 

fc0 _ i_oc 2vr (k + ie) 2 (-k 2 + k 2 + M 2 - ief { ' 

The general power in integral |140| leads to a branch cut along the real axis for all those values for which 
^o — k 2 + M 2 . Including the ie-prescription as indicated in 140, we choose a path in the complex plane around 



the branch cut in the upper half of the plane and find the following solution: 
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v ' J^W+W 2vr k 2 \ - fcg + k 2 + M 2 |/3 

r (l - /?) r (\ + /?) 

= -2 i sin (fin) 3 V . 1 (141) 

2vrl(k 2 + M 2 )5+' 3 

The remaining Euclidean integral can then be performed easily. In the case of the crossed ladder diagram 
Mcl we find in this manner again a divergence which is hidden in Feynman parameters. This can be handled by 
splitting into f(u) and f(l) terms as above. For 



d n k ^ e 7r(/c 2 , m 2 " 1 



we find 



(2tt)» (fco + ie) 2 (A; 2 - A 2 + ie)((k + q) 2 - A 2 + ie) 



e \ r(-l-f)r(l + e)r(2 + f)r ? e yi /-i x(l - x) « 



(142) 



/(«) = 16 sin -7T — ^- \, iLJ ~ \ dx / dv ( 14 3) 

(47r) 4 7rm 2 Jo Jo a 2+ % (A+^u) 



tt^ i« • e \ T (-1 - f ) T (1 + e) T (2 + §) r?- /-i /-i x(l-s) 

/(l) = 16 sm -7T — ^y- ,, — — / dx I dv t— (144) 

V2 7 (47r)%m 2 7o io (ztv(l-v^ + K\ 

The vertex correction graph and the integral occurring in the onle-loop verex correction term M i vc can 

be calculated analogously. Here we have 

f d n k ^(k 2 ,m 2 ) 

J (2tt)» (k + ie) 2 {k 2 - A 2 + ie) 1 j 

with the corresponding solutions 



, e \ r(-f)r(e)r(l + f)r ? e /-i x(l-x) . . 

/(tt) = -16 sm -7T v 27 , /. \ — 2 -— L - / dx H nr (146) 

\2 / (4vr) 4 7rg 2 7 a i+f f i _ u + 4^ 



e \r(-f)r( f )r(i + f)f 



m 



2' 



'<"> - i-l?)^w^^J (147 > 
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